Cooperative order and excitation spectra in the bicomponent spin networks 
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A ferrimagnetic spin model composed of S = i spin-dimers and S = | spin-chains is studied 
by combining the bond-operator representation (for 5 = ^ spin-dimers) and Holstein-Primakoff 
transformation (for S — ^ spins). A finite interaction Jdf between the spin-dimer and the spin 
chain makes the spin chains ordered antiferromagnetically and the spin dimers polarized. The 
effective interaction between the spin chains, mediated by the spin dimers, is calculated up to the 
third order. The staggered magnetization in the spin dimer is shown proportional to Jdf- It presents 
an effective staggered field reacting on the spin chains. The degeneracy of the triplons is lifted due 
to the chain magnetization and a mode with longitudinal polarization is identified. Due to the 
triplon-magnon interaction, the hybridized triplon-like excitations show different behaviors near the 
vanishing Jdf- On the other hand, the hybridized magnon-like excitations open a gap ~ Jdf- 
These results consist well with the experiments on Cu2Fe2Ge40i3- 

PACS numbers: 75-10-Jm, 75.50.-y, 05.30.Jp 

I. INTRODUCTION 

Quantum magnetism has received considerable attention from both theoretical and experimental points of view in 
the past decades. Some low-dimensional magnets, for example, antiferromagnetic spin chains with half odd integer 
spins are gapless and have a disordered ground state; while some others, such as antiferromagnetic spin chains 
with integer spins, spin ladders and dimerized spin chains, are gapped and disordered. In two dimensional cases, 
Hcisenberg antiferromagnets in the square lattice are gapless and have ordered ground states at zero temperature. 
Besides the dimensional effect, various frustration and anisotropy cause novel and complex phenomena. Recently, 
bicomponent systems combining two different spin frameworks have been realized experimentally in R2BaNi03^'^, 
Cu2Fe2Ge40i3^i^i^i^, Cu2CdB206'' and Cu3Mo20g^. Among them, Cu2Fe2Ge40i3, which incorporates intercalated 
Cu^"*" spin dimers [S — i) and Fe'^+ spin chains [S = |), is most extensively studied. Below Tjv = 39 K, a cooperative 
order was observed by the measurements of susceptibility and heat capacity^. At T = 1.5 K, the estimated mcu = 0.38 
fiB and mpe = 3.62 fiB- Compared to the classical expectation value of 1 ^b, "t-Cu is drastically suppressed and keeps 
proportional to mpc at all temperatures. By detailed inelastic neutron scattering study^'-.,6^ Cu2Fe2Ge40i3 was found 
to exhibit two types of spin excitations with separate energy scales. Although the Fe'^~'"-centered low-energy spin 
excitations can be well interpreted by the spin wave theory and a small gap of about 1 meV is estimated, the Cu^'*'- 
centered high-energy part is less understood. Masuda et. al^ guessed the presence of a triplet mode, which should 
have longitudinal polarization and be totally incompatible with conventional spin wave theory. In this paper, we 
combine bond operator representation^ and Holstein-Primakoff transformation— to study this bicomponent system. 

Since the exchange coupling along the c-axis is pretty weak, the real geometry of Cu2Fe2Ge40i3 could be simplified 
as a two-dimensional topologically equivalent model^ (Figd]). The model Hamiltonian consists of the Cu-Cu, the 
Fe-Fe and the Cu-Fe interactions: 

with 



Hqu-Fc 

where Jd and Jdf denote the Cu-Cu and Cu-Fe interactions respectively, and Ja(b) the Fe-Fe exchange constant along 
a(b) direction. By neutron inelastic scattering and neutron diffraction^, Jd was found dominantly larger than other 
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FIG. 1: The geometry of the topologically equivalent spin lattice of Cu2Fe2Ge40i3. The Cu-Cu interaction Jd is dominantly 
larger than the Fe-Fe interaction Ja{b) SLiid the Cu-Fe interaction Jdf- The weak interchain interaction along c-direction is 
neglected for simplicity. 




FIG. 2: The effective four-spin model. 



interactions and Jb is very close to J^. The Fe^+ has spin 5* = | and the Cu^+ has spin T = ^. Single- ion anisotropy 
and external magnetic field are not considered here. With a vanishing Jdf, the model reduces to the independent 
one-dimensional Heisenberg antiferromagnetic chains with S = ^ and isolated spin dimers with T = ^. With a finite 
Jdf, an effective interaction between the Fe^+ chains, mediated by the Cu^+ dimers, will make the Fe'^+ spins stay 
in the Neel state. (The small Fe-Fe interactions along c-direction reinforce this order). On the other hand, the local 
magnetization of Fe'^^ spins will tend to excite the Cu^"*" dimers from the singlet to the triplet and the Cu^+ spins 
will then show local magnetizations. In the following, we study the effective interaction between the Fe'^+ chains 
mediated by the Cu^+ dimers perturbatively and the Cu^+ magnetization induced by the local Fe'^^ spins (section II). 
In section III, we study the excitation spectra of the mixed spin system. A summary is given in section IV. 



II. INDIRECT FE-FE MAGNETIC INTERACTION AND CIJ'^+ MAGNETIZATION POLARIZED BY 

THE FE-'+ SPIN 

As mentioned before, Jd is much larger than Ja(b) and Jdf- Considering the low energy excitations only and 
integrating out the degrees of freedom of the spin dimers, we could obtain an effective Fe-Fe exchange Jeff mediated 
by isolated spin-i dimers^^. In the model IIamiltonian(eq. (1)), J^cu-Fc could be treated as a small amount and i/po 
will be neglected for the moment since no Cu^+ spin operators are involved in it. As shown in Fig. 2, the model then 
reduces to a four-spin system 



n = Ha + H' 
with Ho = JdTi • fa and H' = Jdf(T'i • S'l + fs • S2). 



(3) 
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Hq has a singlet ground state \s > with a ground state energy Eq = — |Jd- Its excited states are triplets 
I^Q, > with a — 0,±1, corresponding to the total z component Ti^ + = 0,±1. The eigenenergy of the ex- 
cited states is Ei^a = jJd- The first order perturbation of iJcu-Fc is Ti-i —< s|iJcu-Fc|s >= 0. The 2nd- 

order perturbation = Ea <"l^'^— ^°l^°><*°l-^'""^'"''l''> = ^5-^.5-2. The 3rd-order perturbation is = 

X;„,/3 <"l^^"-'-° l*° ><^°j^^j--^^j[*^>g ^ T^Si ■ S2. Thus, up to the third order, the effective exchange 

coupling mediated by the spin-f/2 dimer is obtained as 

Ti-cS = JcffSi ■ S2 (4) 

with Jeff = §^(f + ^). 

With the effective interaction between the Fe^+ spin chains considered, the Fe'^+ spins will order antiferromagneti- 
cally. We use a molecular field approximation to study the effects of the staggered Fe'^+ magnetization on the Cu^+ 
spin dimers. It is relevant at a higher energy range (comparable to the dimer gap), where the dynamics of the system 
is dominated by the dimers and the chain freedom can be effectively integrated out. The four spin Hamiltonian 
becomes 

H = Jufi ■ f2 + Jdf"1fc(T'iz -T2z), (5) 
where mpc — < Sz > is the staggered magnetization of Fe'^+ spin. This approximate Hamiltonian can be exactly 
diagonahzed with the bases \s > and \ta > (a = 0,±f). The eigenvalues are eg = —jJd — ' ^"^"''™^" 



ei = —jJb + l-^Dy 1 + ' dj"'f° ^ and 62 = 63 = jJb- The corresponding eigenstates are \^pi >= a^js > +bi\to > 

{i = 0,1) with a, = _^'T7:r2 2 and h = /, J'^I'L^ 2 ^ 1^3 l^i > and |^4 >= \t-i >. From the 

ground state \ipo >, we find that \to > is partly excited and the Cu^"*" spins now get a finite staggered magnetization: 
mcu =< T2Z >= — < Tiz >=< '0o|T2z|'0o >= oo&o- When mpc or is small, we have rncu ~ ^^wpc, which agree 
well with the experiment^. 

III. EXCITATIONS IN THE BICOMPONENT SYSTEM 

In the Fe-Cu spin system, if we neglect the interaction between the Cu^+ and Fe"^+ spins, the Fe'^+ spins will be 
in the Neel state (considering the antiferromagnetic interaction along c-direction) and its excitations are spin-waves 
(magnons); while the Cu^+ dimers are in the singlet states and the excitations are triplets (triplons). When the Cu-Fe 
interactions are switched on, the Neel order in the Fe"^^ spin network is reinforced and a staggered magnetization in 
the Cu^"*" spins are induced. The magnon and triplon interact with each other and the hybridized excitations show 
interesting behavior. In the following, we study these excitations. 

For Cu^+ spin dimers, we use the bond operator representation. With the definition of \s >= s^|0 >, |to >= ^^|0 >, 
>= (it|0 >, \ti itt|0 > and s''s + tH + dU + u^u — 1, the Cu^+ spin operators are expressed as^: 

T+ = -^[-{uls,-sld,)+tld,+uUa]: 

-^[{uls,- sld,)+tld„ + ult„]\ 
T: = ^[{siu + tls„) + ulu, - did,]; 
TUi = l[-{sit, + tis,) + uiu,^did,]; (6) 

where, ct = 1,3, corresponding to the pairs of {Ti,T2) and (TsyT^), respectively. 

We apply Holstein-Primakoff transformation to Fe^+ spinsi^. Supposing 6*1 and ^3 point up and 6*2 and ^4 point 
down, we have 



5*+ = \/2S - ajaiai, 



S, = a]\/2S - a\ai, 



St = 5-a|a„ z = 1,3; (7) 
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and 



S; = b%-S, J = 2,4. (8) 

Before we substitute these representatfons into Hamiltonian (1), some approximations have to be made. As found 
in section II, the Cu^+ dimers are in a mixed state of |s > and \to >. We introduce another four operators asi^^ 



Xrr \ ( COS0 —sm9 



Y„ \ sm9 cos6 \ t, 



Sn 



(9) 



and 

(10) 



\ _( cos6 sinO \ f Xa 
t„ ) ^ \ -sinfl cos6' J \Y„ 

with cr = 1, 3 and suppose bosons are condensed with < Xa >—< X\ >= 1, which means a fong range order in 
the Cu^+ spin network. 9 will be determined variationally or by canceling the single-operator terms. 

For the Holstein-PrimakofF transformation, we empfoy the usual linear approximation: Sf = y/2Sai, S[' — \/2Sal, 
5f = S* — ajfli, i — 1,3 and S*^ — \/2Sbj, SJ^ — V^Sbj, Sj = —S + b^jbj, j — 2, 4. Substituting these transformations 
into Hamiltonian (1), we get 

H = ^[(Jdf5cos20 - i JDsin20)(y/- + + + Fg,-) + (JDF5sin20 + JDCOs20)(y/-Yi,- + ygV^ar) 

r 

+ - Jd(1 + cos26){u\pUip + d\pdip + ul^usp + dlpd^p) + (2 JaS* + - JDFsin20)(aJ-aiF + 52^&2r + Og-osf^ + blpb^p) 



2S TT 

— JBFSm{9 + -){alpulp+ a^^-Wg- + 52rdlf + fo4rC?3f + h.c.) 
2*S 7T 

^ JDFsin(6' - -)(alpdif + al^dsf + ulpb2f + ulpb4f + h.c.) 

+JaS{a3fb2f+ h.c.) + ^JiiS{aifb^p_g + a'^^_^^^blp+ h.c.) + ^Ji,S{aifb2f+s + a3fb4f-a + al^^^bl^ + al^^^blp + h.c.)] 

-2A^JDF'5'sin26' - + cos26') Jd - ANJ^S"^ + 3 - operator terms + 4 - operator terms. (11) 

The 3-operator and 4-operator terms are omitted here for simplicity. By letting the coefficient of the single-operator 
terms be zero, we get 

9 = -arctan — 2^— . (12) 
2 Jd ^ ' 

Correspondingly, at zero temperature, we get the staggered magnetization at the Cu site as 

mcu = I < Ti, > I - -sin20 = , (13) 
2 ^i + (2iQ£5)2 

When Jdf << Jd, 9 ~ -'^J'^ and tocu ~ ^^j^- This result is consistent with mcu ~ ^^wpe, which we get from the 
molecular field approximation in section II. 

After Fourier transformation, the Hamiltonian can be written as 

H = Y.ujl{Yl,y4,k + Yl^Y3^u) + ^lHk^k + C, 

k 

where, 

LoX ~ JdySs\d29 + Jucos29, 
C = -2NJuFSsm29- N{^+cos29)Jb-4:NJ^S^, 



*fc = ia{,k^b2-k,aU,b4-k,d\^^^,ui-k,dU,U3-k), 



5 



H,, = 



/ 


a 


Se 







— 7C0SJ^ 


7sini/ 





\ 




Jg'ifcy 


a 


5 





7sini/ 


— 7C0SJ^ 













5 


a 


^gjfca 








— 7Cosi^ 


7sinj/ 











a 








^smu 


— 7Cosi^ 



























7sin:/ 


— ^COBV 











d 








V 










7sinj/ 








d 











7siiw 


— 7C0Si^ 











d } 



(14) 



with 



2Ja5+ ijDFsin26l, 
ijD(l + cos26'), 
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— JaS, 















4 



(15) 

Compared with the conventional spin wave theory, the magnon gets an additional energy of ^ JDFsin20 = Jdf^Cu, 
which comes from the polarization of the Cu^^ spin dimer. As a result, the magnon excitations show a tiny gap. We 
will discuss it later. The Hamiltonian can be diagonalized as 

k 1=1,8 k 

where, a constant has been neglected and the expressions for ujik are given in Appendix A. 

The triplon excitations Yo-fc {a = 1, 3) do not interact with the magnons. They are dispersionless. The spectrum 

Jd[1 + 2(^^7^)^]/y/l + 4(^^aE£)2 ^ j^^i 4^ 2(^^Q^)'*]. It increases with increasing Jdf and returns to the 

singlet-triplet gap when Jdf = 0. With the Cu^+ spin operator (eq. (6)) expressed by the and Ya (eq. (10)) and 
making mean- field approximation of < >= 1, we find proportional to Y and Y^ and (T'^iv)^ proportional 
to u{u^) and d{d^). Thus the operators Ycrk = sin^Scrfc + cosOtcrk {a — 1, 3) describe the longitudinal fluctuations and 
we believe they are the modes having longitudinal polarization and totally incompatible with conventional spin wave 
theory, as guessed by Masuda et. al^. It is emphasized that these modes cannot be obtained by spin wave theory or 
effective spin wave theory^. 

Among the eight aik excitations, four branches are magnon-like and another four are triplon-like. We denote 
them by Aik, A^k, B2k, B^k and Dik, -Ds/t, Uik, U^k respectively. The corresponding spectra are represented by 
w^^, U!^^ , Lo^^ and ^k^ ^ '^k^ ^ '^k^ ■ They are two-fold degenerate. At Jdf = 0, they return to the pure 

magnon spectra uj^^^ — Ja.S\J 2(1 ^ cos-^) and pure triplon spectra w^^^ = Jd. The magnons are gapless and the 

triplons are dispersionless. At a finite Jdf, the magnons become gapful and the triplons become mobile due to the 
magnon-triplon interaction. 

Experimentally, Jd and Ja were determined as Jd = 24 meV and Ja = 1.6 meV. The reported Jdf is discrepant 
such as 2.4 meV, 0.9 meV or 2.0 meV. In Fig. 3 and 4, we show the numerically calculated magnon-like and triplon-like 
excitation spectra, respectively, with Jd = 24 meV, Ja — 1.6 meV and Jdf = 1 meV. The calculated magnon gap 
is 1.27 meV. Experimentally, a small empirical anisotropy gap is estimated as 1 meV (spin wave theory-) or 2.02 
meV (effective spin wave theory^). We find that the gap largely depends on Jdf- The triplon-like excitations have a 
dispersion with small amplitude. 

Regarding ^^j^ as a small amount, we could get analytical expressions for the spectra. The details are given in the 

Appendix. The magnon-like excitation has a gap — a;"^''(0,7r) — As — lo^"' {211,11) w Jd\/ j^^\p\ (o" = 1,3, 



a' = 2,4) with ry = ^nd p = —^^f^\J f ■ It should be pointed out that the magnon-like spectra are proportional 
to {^^j^Y away from the momentum of (0,7r) or (27r,7r) (see eq. (jAQp ). The gap of the triplon-like excitations is 



A[/ = w^'' (27r, 7r) = Ad ~ (jJ'°"(0, tt) « Jd(1 — jzh^P^) (^ = li3). When rj « 1, the gaps can be further simphfied 
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FIG. 3: Spectra of the magnon-like excitations '^k^ y y ^k"' with Jd ~ 24 meV , Ja = 1.6 meV and Jdf ~ 1 meV. 

Compared to the case of Jdf = 0, a smaU gap opens. 
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FIG. 4: Spectra of the triplon-hke excitations uj^ , '^'j^" and uj^" {a = 1,3) with Jd — 24 meV , Ja — 1.6 meV and Jdf = 1 
meV. The dispersionless oj^ is longitudinally polarized, as mentioned by Masuda et. al.— Compared to the case of Jdf = 0, 
oj^" and uj'^" {a — 1,3) show weak dispersions. 



as Aa « JoF^'y ^ « Jdf and Au « Jd[1 - j^i'^^^)'^] i"^)^- These relations are numerically shown in Fig. 
5. Reminding Ay cx (^^j^)"*, we find that different branches have different dependences on The band width of 
the triplon-like excitations is Wu = Wd ~ Jdjz^P^ ~ Ja( "^jp"^ )^- As a counterpart, an effective interaction between 
the Cu^+ dimers, mediated by the Fe'^"'" chains, can be estimated as J^g '--^ Ja( '^°^'^ )^. 

Experimentally^iS., the magnetic moment in Cu2Fe2Ge40i3 was found to be nearly confined in the a-c plane. 
Dzyaloshinskii-Moriya (DM) interaction and other insignificant anisotropy effects may exist. For simplicity, we neglect 
these anisotropy effects and the weak interaction along c-direction in our mean-field theory. Some remarks on the 
Goldstone theorem have to be made here. We start from an isotropic spin Hamiltonian and a gapless mode is expected 
in the ordered state because of the rotational symmetry. The small but finite magnon gap seems unwelcome. We 
interpret it by a two-step process. Similar to the mixed-spin antiferromagnets R2BaNiOg^, there is a separation of 
energy scales of magnetic excitations in Cu2Fe2Ge40i3. The Fe'^^ centered magnons with low frequencies and the 
Cu^"*" dimer centered triplons with high frequencies have different dynamical behaviors and different timescales. At 
the mean-field level, we first consider the effective Fe-Fe interaction Jcg and the induced long range order in the 
Fe'^+ sublattices. Goldstone modes appear at this stage. We then study the polarization of the Cu^+ dimers and its 
reaction on the Fe'^"'" sublattices. The hybridized magnons get a small gap due to the staggered field presented by 
the polarized Cu^"*" dimers. This phenomenon has been studied earlier in spin-i and spin-1 antiferromagnetic spin 
chains experimentallji^ii^ and theoreticalljsi^. It is shown that the presence of a staggered field will make the spin-i 
antiferromagnetic chain gapful and split the Haldane triplet into two branches. In our studies, we further consider 
the fluctuations of the staggered field and their hybridization with the magnons. It may be interesting to consider the 
DM interaction or other anisotropic effects, which break the rotational symmetry and lead a gap naturally. We argue 
that the polarized Cu^+ dimers will make important contributions to the gap, even dominant if the anisotropies are 
tiny. We note that the estimated gap is close to the estimated Jdf from either the spin wave theory^ or the effective 
spin wave theory*. 

In the case of Jdf < 0, the Cu^"*" spin T„ will be parallel to its nearest Fe'^+ spin Sn [n = 1, 2, 3, 4). Our theory is 
still applicable. The angle becomes 6' = — ^arctan^i^^^^. In Fig. 5, we also show the dependence of the energy 
gaps on in the range of Jdf < 0. 
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FIG. 5: The magnon-like excitations have a gap Aa proportional to Jdf (squares), while the triplon-like excitation shows Jj^p 
{Ay, triangles) or Jpp dependences {Au, circles). Jd ~ 24 meV and Ja = 1.6 meV. When Jdf < 0, the Cu^^ spin is parallel 
to its nearest Fe^"^ spin; while when Jdf > 0, they are antiparallel. 



IV. SUMMARY 



In summary, we combine the bond-operator representation and Holstein-Primakoff transformation to study the 
mixed spin lattice model made up of Cu^"*" {S = i) spin-dimers and Fe'^^ {S = |) spin chains. A finite interaction 
Jdf between the spin-dimer and the spin chain introduces a cooperative order. The effective interaction between 
the spin chains, mediated by the spin dimer, is calculated up to the third order. It makes the Fe'^+ spins ordered 
antiferromagnetically. Meanwhile, the local magnetization of Fe^+ spins polarizes the Cu^+ spin dimers and the 
staggered Cu^"*" magnetization is shown proportional to Jdf as well as the local Fe^+ magnetization. It effectively 
presents a staggered field reacting on the Fe^+ spin chains. Considering the triplon-magnon interaction, the excita- 
tion spectra are especially investigated. The degeneracy of the triplons is lifted and the mode having longitudinal 
polarization and totally incompatible with conventional spin wave theory, as pointed out by Masuda et. al, is identi- 
fied. It shows ( '^jp'^ )^ dependence near vanishing Jdf- The hybridized triplon-like excitations get a weak dispersion 
and show a different (^^)^ dependence. At the same time, the hybridized magnon-like excitations open a gap 

— Ab ~ JufS 'fj^^ ~ Jdf at the momentum of (0, tt) or (27r, tt). Away from this momentum, the magnon-like 
spectra show (^^)^ relation, instead. The experiments on Cu2Fe2Ge40i3 are interpreted. 
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APPENDIX A: EXCITATION SPECTRA 

From the matrix Hk (eq. (|14p ). we get the equation to determine the excitation spectra: 

to^ + Uuj^ + Vlj^ + Wu? + Z = 0, (Al) 

where, 
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U = -2(0^ + (f -26"^ + 2j'^cosip), 

y = + 4a^d^ + d^- iad-y"^ + 87^ - 4a'^5^ - Sd^S"^ + 2S^ - 2S^cosk^ + 4'f{o? + d^ - 2S'^)cosip 
+37''cos^93 — Adj'^ScoskySmip — Sj'^sm^ip, 

W = -2a^d' + Aa^'dj^ + Ad^S"" + 27^52 + 40^72(^2 _ 2^2) _ a'{2d^ + 7^ - 8^252) _ ^2(^4 ^ 4^4^ 
+Ad5^coskx{d5 + ^"^coskyshup) + 7^(-(7^(a2 + rf^ - 252)cos(2^)) - 7''cos(3<^) 
+Ad5{a^ + c/2 - 6'^)coskyS\rnp + cosip{-Aa'^d^ + Sad-y"^ - 87^ + Sd'^S'^ + Sd-y'^ScoskySijup)) 

Z = aW" - 40^^37' + ^ - + 2d"<5' + 02(5^27^ - AdH^) + 2a(-d76 + 4^^7252) 

8 2 

-(i2<52cosfc^(74 + 2^2 52 _ ^4pos(2(^)) + ^7^(4(202^2 - 40^72 + 7^ - 5'^)coa{2>^) + 7''cos(4(^)) 

8 

-2d72(5cosA;y(2o2d2 - 40^72 + 7^ - 2d^5'^ + 2d'^5'^cosk^ + 7''cos(2(^))sin(^ + d^'y^5'^{l + 2cos{2ky))sm^ip, 

with ip = 2(6*+ f). 
Its roots are 



= -\{U + e^/K + m-AV) + <j^^{U + e^/K + U^-AVy - + e-j£L^==), (A2) 
where e = ±1, u = ±1 and 

^=^(4-/(If^)"-(-l-/(l)MIf)''^ (-) 

with p = — j^ — Z and q = — — ^ + If | Jdf| J^d, the triplon-hke excitation spectra can 

be expanded as 

Ji' = JD(l + 5i(fcy 

wf" = JD(l+5-i(fcy), <7=1,3, (A4) 

where, p = —^f^\f^ and 

7?(^ - 4772 + 8 - cosA;x(8?72 - 4 + rj^cosky) + 2e(?72 + (772 - l)cosfcx)-v/3Tcos(2fc^j 
^^^^^ = -r74 + 8r?2_8 + r7Wj, ^^^^ 

with 77 = ^^p- The energy gap of the triplon-Uke excitations 

A[, = Az5 « Jd(1 - Y^^-^')' 

locating at the momentum of (27r, tt) (for the branch of u)^" , cr = 1, 3) or (0, tt) (for the branch of w^' , f = 1, 3). The 
band width of the triplon-hke excitations is 

Wu = Wd^ Ju t^p" ■ (A7) 
I — r]-' 

The energy gap of the magnon-like excitations occurs at the momentum of (0, tt) (for the branch of co^" , cr = 1, 3) 
or (27r, tt) (for the branch of co^" ,(7 = 2,4) and shows linear dependence on Jdf: 



Aa = Ab = Jd^/ jIpI « I Jdf|<S'\/ ^— I » I^dfI- 



(A8) 
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Away from that momentum, the magnon-hke excitation spectra can be expanded as 




a 



= 1,3. 



(A9) 



Interestingly, they exhibit quadratic dependence on ^W- . The band width of the magnon-hke excitations is 
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